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We calculate theoretically the Coulomb drag resistivity for two graphene monolayers spatially 
separated by a distance "d". We show that the frictional drag induced by inter-layer electron- 
electron interaction goes asymptotically as T^/n"^ and T'^\n{n)/n in the high-density (kpd ^ 1) 
and low-density (kpd <C 1) limits, respectively. 
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Frictional drag measurements of transresistivity in 
double layer systems have led to significant advances in 
our understanding of density and temperature depen- 
dence of electron-electron interactions in 2D systems [l| . 
Recent interest has focused on the role of electron in- 
teraction effects on the graphene drag resistivity, which 
should vary in a systematic manner as a function of elec- 
tron density (n), layer separation (d), and temperature 
(T). In particular, a recent experiment of Coulomb drag 
in double layer graphene by Kim et al. 2] is particu- 
larly interesting. In view of the considerable fundamental 
significance of the issues raised by the experimental ob- 
servations, we present in this paper a careful theoretical 
calculation of frictional drag Pd{T) in a 2D [3*1 graphene 
within the canonical many-body Fermi liquid theory. The 
current work is a generalization of the earlier theoretical 
work on graphene drag by Tse et al. 0] ■ 



screened electron-electron interaction is given by 



We start by writing down the theoretical formula for 
in the many-body Fermi liquid RPA-Boltzmann 
theory approximation widely used in the literature. The 
double layer frictional drag in a many body-Fermi liq- 
uid diagrammatic perturbation theory with dynamically 
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where Fi^2{q,oj) = |u^^(g,tj)|Imrii,22(q,a;), with 
^12 ~ ^12/ ^(It^) dynamically screened interlayer 

Coulomb interaction between layers, and T{q,uj) is the 
2D graphene non- linear susceptibility [4]. The interlayer 
Coulomb interaction is given by Vi2iq) — v'^(q)e^'^'^ with 
the intralayer Coulomb potential ^Sj] Vc{q) = / nq 
where k is the background dielectric constant. (We con- 
sider the so-called balanced situation here with the same 
carrier density n in both layers.) Note that the dielectric 
function e[q,ui) entering Eq. ([1]) is given by '3] 

= [1 - «ll('7)nil(g,w)] [1 - V22{q)^22{q,i^)] 

- ^12(9)1^21 (9)1111(9,^)1122(9,0;), (2) 

where Va^q) — vdq) and 11^^ is the intralayer graphene 
polarizability Even tough there are analytic expres- 
sions for graphene polarizability at T = 0,[5| the finite 
temperature versions of the polarizability have not been 
calculated analytically. The full expression of finite tem- 
perature polaizability is necessary to understand more 
precisely the temperature dependent drag including the 
plasmon enhancement effects [6| . Here we provide the effi- 
cient way for calculating the finite temperature graphene 
polarizability by generalizing our earlier work j5|] 
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where H = H/Z?o (-Do = 2kF/TrfiVp is the density of 
states of graphene at Fermi energy, and kp and vp are 
the Fermi wave vector and Fermi velocity of graphene), 
q = q/kp, OJ = u:/Ep, f{q) - [g-^^^"^) - g{q) = 

[e(e9+M) + hvpq, and a± = up- - (f' ± lujk. 

In /(g) and (7(9), p is the finite temperature chemical 



potential which must be calculated self-consistently to 
conserve the total electron density. 

With assumptions of a large inter-layer separation 
(kpd ^ 1, or qTpd 3> 1, with qxp being the Thomas 
Fermi (TF) screening wave vector) and the random phase 
approximation (RPA) in which Ha is replaced by its value 
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FIG. 1. The temperature dependence of Coulomb drag for 
(a) d = 50A and (b) d = 200A . The black (red) lines are the 
results calculated with equal electron densities, m = n2 = 



XW^cm,-'- (m = n2 = 2 X 10' 
drag resistivity pd(T)/T'^. 



Insets show the scaled 
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for the non-interacting electrons, we have for the identi- 
cal layers at high density and low temperature. 
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where qrp = '^rgkp is the TF wave vector with the 
graphene fine structure constant = e'^/nhvp and C 
is the Riemann zeta function. This result shows that 
Poin) oc n~^, but this result, which was obtained in 
ref. [4], applies only for high density and large separa- 
tion limits (or weak interlayer correlation, kpd ^1). At 
low densities (or strong interlayer correlation, kpd 1) 
the exponent in the density dependent drag differs from 
-3 as shown in Fig. [2] where we directly numerically cal- 
culated graphene drag using Eqs. (Pl-dll). Eq. |4] shows 
Pd{T) oc T^. For large layer separation (i.e. kpd ^ 1) 
the back-scattering q 2kp is suppressed due to the ex- 
ponential dependence of the interlayer Coulomb interac- 
tion vi2{q) oc exp{—qd)/q as well as the graphene chiral 
property. In this case the drag is dominated by small 
angle scattering and one expects po oc T'^/{n^d'^). 

For the strong interlayer correlation {kpd 1) in the 
low-density or small-separation limit, we have the follow- 
ing asymptotic behavior: 
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Thus, we have Pd{T) ~ again, and poin) ~ ln(n)/n 
with very weak logarithmic d-dependence. 

In Fig. [1] we show the calculated Coulomb drag as a 
function of temperature for two different densities n = I, 
2 X lO^^CTO"^, layer separation (a) d = 50 A (h) d — 
200 A. The overall temperature dependence of drag is 
close to the quadratic behavior. But we find a small 
corrections, especially at low values of kpd. In regular 2D 
systems there is a ln(T) corrections to the dependence 
of the drag. However, due to the suppression of the back- 
scattering in graphene such logarithmic correction does 



FIG. 2. The density dependent Coulomb drag for different 
layer separations d = 5, 10, 20 nm and T — 200K. Dashed 
line indicates po ~ behavior. 



not show up in our numerical results except perhaps at 
extremely low temperatures. 

In Fig.[2]the density dependent Coulomb drag is shown 
for different layer separations. Our calculated Coulomb 
drag resistivity follow a n" dependence with a < — 2 
at low carrier densities (or, kpd < 1), but as the den- 
sity increases the exponent (a) decrease. Based on our 
calculation we believe that the experimental departure 
from the behavior reported in Ref. is essentially a 
manifestation of the fact that the asymptotic regime 
is hard to reach in low density electron systems where 
kpd 3> 1 limit simply cannot be accessed. We predict a 
weak ln(n)/n density dependence in the low-density or 
small separation limit. 

In conclusion, wc study the frictional drag between two 
spatially separated graphene layers within a many-body 
Fermi liquid theory. We find that the temperature depen- 
dent drag mostly shows a quadratic behavior regardless 
of the layer separation, but the density dependence varies 
from ln(n)/n for kpd <C 1 to for kpd ^ 1. But most 
currently available double layer graphene samples belong 
to kpd ^ 1, so the density dependence dose not have any 
universal power law behavior. We also find that due to 
the suppression of the q — 2kp back-scattering there is 
no ln(T) correction in the drag resistivity in graphene. 
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